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ABSTRACT 

With the aim of constructing accurate 2D maps of the stellar mass distribution in nearby galaxies 
from S 4 G 3.6 and 4.5 fim images, we report on the separation of the light from old stars from the 
emission contributed by contaminants. Results for a small sample of six disk galaxies (NGC 1566, NGC 
2976, NGC 3031, NGC 3184, NGC 4321, and NGC 5194) with a range of morphological properties, 
dust contents and star formation histories are presented to demonstrate our approach. To isolate 
the old stellar light from contaminant emission (e.g. hot dust and the 3.3 fim PAH feature) in the 
IRAC 3.6 and 4.5 /im bands we use an Independent Component Analysis (ICA) technique designed 
to separate statistically independent source distributions, maximizing the distinction in the [3.6]-[4.5] 
colors of the sources. The technique also removes emission from evolved red objects with a low mass- 
to-light ratio, such as asymptotic giant branch (AGB) and red supcrgiant (RSG) stars, revealing 
maps of the underlying old distribution of light with [3.6]-[4.5] colors consistent with the colors of 
K and M giants. The contaminants are studied by comparison to the non-stellar emission imaged 
at 8/xm, which is dominated by the broad PAH feature. Using the measured 3.6/im/8/im ratio to 
select individual contaminants, we find that hot dust and PAH together contribute between ~5-15% 
to the integrated light at 3.6 /zro, while light from regions dominated by intermediate-age (AGB and 
RSG) stars accounts for only 1-5%. Locally, however, the contribution from either contaminant can 
reach much higher levels; dust contributes on average 22% to the emission in star-forming regions 
throughout the sample, while intermediate age-stars contribute upwards of 50% in localized knots. 
The removal of these contaminants with ICA leaves maps of the old stellar disk that retain a high 
degree of structural information and are ideally suited for tracing stellar mass, as will be the focus in 
a companion paper. 
Subject headings: 



1. INTRODUCTION 

Accurate maps of the stellar mass distribution in 
nearby galaxies are essential for charting the dynamical 
influence of stellar structures over time. They provide 
a key perspective on the evolution in the amount and 
distribution of baryons from the high redshift to the lo- 
cal Universe, both as the result of internal, secular pro- 
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cesses (e.g. Kormendy fc Kennicutt| |2004 Elmegreen 
et aL| |20"07l |Zhang k Buta||2007| jHaan et al.||2TO| 



Foyle et alpOlOl |Elmegreen et al.||2009| |Elmegreen fc 
Elmegreen 2010), and given 'bottom up' assembly and 



external (heirarchical accretion ) events within the con 
text of CDM cosmo l ogy ( e.g., [Balcells et al. 2003 and 
|Courteau et al. 1996; Fl orido et al.|2001[ and see Navarro 
Ife White||1994t |Somerville||20021 |Governato et al.||20071 
|Guq fc White||2008| ). 

Yet, stellar mass estimation as a field is still maturing. 
NIR bands are often exploited as optimal windows on the 
old stars that dominate the baryonic mass in galaxies, 
with images at these wavelen gths serving as good rel- 
ative proxies for stellar mas s ([Elmegreen & Elmegreen 
1984| Rix fc Zaritsky||1995| |Grosb0l||1993p. But esti 



mates incorporating dependencies on star formation his- 
tory, me tallicity, and especially the age of the stellar pop- 
ulation ( |Rix fc Rieke|1993||Rhoads|1998[) , as well as dust 
extinction and reddening, have only relatively recently 
emerged, following the demonstrated correlation between 
optical colors and the m ass-to-light ratio in stellar popu 
lation sy nthesis models (|Bell & de Jong||2001 
al ~ 

Spitzer Survey o f Stellar Structure in Galaxies 



latio n sy : 
aT1[2^09l 
^TEe~Si 



Zibetti et 



(S G; Sheth et al. 2010) provides an un-paralleled in- 
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ventory of stellar mass and structure in nearby galax- 
ies, with deep imaging at 3.6 fim and 4.5 /im (reach- 
ing fJhi.&iim{AB) (1(t)~27 mag arcsec -2 ) where the light 
mainly traces the oldest stars. The sample of S 4 G galax- 
ies therefore stands to significantly contribute to our un- 
derstanding of the evolution of stellar mass and struc- 
tural properties. 

For the purposes of accurately mapping the radial dis- 
tributions of mass in galaxies, and in particular later 
Hubble types, the situation remains less than ideal. Al- 
though M and K giants should in fact be the foremost 
contributors to the flux in these two shortest IRAC wave- 
bands, by now several other sources alongside the old- 
est stellar light are well-discerned at IRAC wavelengths 
in star-forming galaxies. The finite instrumental band- 
width introduces additional features (i.e. the 3.3 \xm 
PAH emission feature in the 3.6 /xm ban d and the PAH 
continuum at 4.5 fim; Flagey et al.|2006 ), accessory non- 
stellar continuum emission (i.e. from hot dust) appears 
at these wavelengths, and lower-M/L asymptotic giant 
branch (AGB) and red supergiant (RSG) stars can also 
make significant contributions to the emission. 

Depending on physical conditions in the disk, the de- 
gree of active compensation for these sources varies from 
galaxy to galaxy, and different investigators propose 
different compensation schemes. Total mass estimates 
should be only modestly affected by the presence of con- 
taminants, which contribute a few percent to the inte- 
grated flux. But ID and/or 2D representations of either 
the relative or absolute stellar mass estimated from 3.6 
/im images can be greatly influenced by the pronounced 
local contribution made by these sources. These three 
main contaminants can lead to over-estimation in stel- 
lar mass maps and therein complicate measurements of 
structural diagnostics. 

Several schemes for the identification (if not correc- 
tion) of these features in IRAC images have be e n devel- 
oped JKendall et al.|[2008l |Hunter et al.||2006| |Bolatto 
et al.||2007| |de Biok et al.||2008D, but only few with the 



aim of deriving maps of the mass distribution in nearby 
galaxies. Given that PAH and hot dust contaminants 
appear in star-forming (high-density) spiral arm regions, 
distinct components fits to the stellar emission there are 
prone to inaccuracy and so are typically omitted in favor 
of an overall azimuthall y-averaged contribut ion insur 



face brightn ess models ( Kendall et al. 2008| and Leroy 
et al.||2008|). While these methods can often supply a 



contammation-free view of the old stellar disk, it is not 
immediately clear that they can be successfully extended 
to a large data set like S 4 G, which samples throughout 
the Hubble sequence, over a large range of dust contents, 
gas fractions, and star formation histories. Furthermore, 
the auxiliary data employed by these methods (e.g. 8 
/im, UV, optical images) will not be uniformly available 
for the full S 4 G sample of 2,331 galaxies. A correction 
for contaminant emission at 3.6 /im and/or 4.5 /im that 
is self-reliant, and can assure uniformity throughout the 
sample, is therefore highly desirable. 

Corrections for non-stellar emission are especially crit- 
ical for mass estimates, since this emission is not incor- 
porated into SED libraries at the core of now-st andard 
techniqu es for estimating th e stellar M/L (e.g. |Bell fc| 
de JongpOOT] |Zibetti et al.||2009| ). Tight constraints on 



the stellar M/L are possible using optical or NIR col- 
ors that leverage metallicity, age and SFH variations in 
the underlying stellar population. Recently, Zibetti et al. 
(2009) have shown that a second, optical/NIR color can 
significantly boost the pred ictive power of the colo r-M/L 
relations first assembled by Bell & de Jong (2001 1. With 
color maps replacing global colors, plus a prescription to 
account for the effects of dust extinction/reddening, lo- 
cal estimates of the M/L in NIR bands can be obtained 
with high accuracy. The few existing 2D maps of the 
stellar mass distribution derived in this way show great 
promise, and reveal a unique glimpse of the changing 
appearance of stellar structures between light and mass. 
Paired with similar information, S 4 G images free of emis- 
sion from contaminants should be readily convertible into 
high-q uality mass maps, as considered in Paper II of this 
series ( |Meidt et aL||2011[ ). 

To complement s 4 G science, we therefore introduce a 
technique that removes contaminants from S 4 G images 
to reveal an underlying uniformly old distribution of stel- 
lar light, with [3.6]-[4.5] c olors that are cons i stent with 



al. 


2006b 







nique used to pursue these ends (described in the next 
section) in principle allows much more information in in- 
dividual structural components, like spiral arms, to be 
maintained in contaminant-corrected mass maps than is 
supplied by other commonly-employed mass estimation 
methods; models of the 3.6 /im and/or 4.5 fim emission 
(wherein the appearance of contaminants is reduced on 
azimuthal averaging), have so far provided the basis for 



most mass estimates at thes e wavelengths (e.g. de Blok 
et al. 2008 Leroy et al.|2008 ) . Our maps of the old stellar 



ight, with full structural information, are ideal relative 
mass tracers for opt imal measurement of bar strength 
and ellipticity (e.g. Menendez-Delmestre et al. 2011), 
spiral arm-interarm contr asts and relative Fourier am- 
plitudes and phases (e.g. lLaurikainen et al.||2004[ |Salo 



et al.|2010b||Elmegreen et al.|201ip . Likewise, we antici 



pate that these maps can be converted into high-quality 
absolute mass maps relevant for deriving the underly- 
ing stellar p otential, measuring bar and/or spiral arm 
tor ques ( e.g. |Stark|| 19771 |Zaritsky & Lo l 1986 Blocket 



al.| |2004| iHaan et a l. 2009; Salo et al. 2010b) and mass 
flux / gas inflow (e.g. Quillen et al. 1995; Zhang & Buta 



20071, and for constraining the baryonic cont ribut ion to 
obser ved r otation curves (e.g. de Blok et al.[ 2008) 



In § 2.2 we describe our approach (motivated in § 2.1 ), 



with particular emphasis on the use of 3. 6 an d 4.5 /im 
S 4 G data and data products (§§ 2.3 and 2.4). Results 



from the application of ICA to a small sample of galax- 
ies (described in §[3| using S 4 G archiv al images processed 
through the S 4 G reduction pipeline (Sheth et al. 20101 
are presented in § |4j These ICA sources are identified 
as either old stars or contaminants, where the latter is 
defined as all sources of emission apart from the oldest 
stars: PAH, hot dust, and intermediate-age red super- 
giant (RSG) and asymptotic giant branch (AGB) stars. 

In §p]we focus on the contaminant emission and relate 
the intensity and spatial information contained in our 
ICA source maps to archival IRAC channel 4 images, 
which are dominated by the broad 7.7 /im PAH feature. 
This comparison reveals the detection of both the 3.3 /im 
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PAH feature and hot dust in 3.6 [im images with ICA (§ 
5.1 1. It also allows for the identification of a second, non- 
dust contaminant at 3.6 and 4.5 /zm, namely the bright 
but low M/L red obj ects contributing to the emission at 
these wavelengths (§ 5.2). 

The collective removal of these contaminants with ICA 
reveals maps of the underlying old stellar light that 
should be readily convertable into mass maps. In a forth- 
coming paper we will examine our ICA maps of the old 
stellar light in detail, and compare different techniques 
for assigning the mass-to-light ratio at 3.6 [im. Here we 
include a preliminary assessment of the quality of these 
maps (as well as of the contaminant correction, itself) 
by comparing the old stellar emission at 3.6 fim with 
other NIR tracers. Finally in § [7] we summarize our 
findings and describe avenues for future work with the 
ICA-idcntificd sources. 

2. ON THE APPLICATION OF ICA FOR 
CORRECTING S 4 G IMAGES 

2.1. Motivation 

Mass maps for much of the S 4 G sample will neces- 
sarily rely on information from only the two shortest 
IRAC wavebands available during the Spitzer warm mis- 
sion. Given this constraint, the task of removing con- 
taminants from these images might involve some sort of 
scaling and subtraction of the two measurements based 
on an (S 4 G-defined) empirical recipe. 19 On the other 
hand, the constraints also suggest an approach that is 
similarly-premised and yet less susceptible to empiri- 
cal interpolation-uncertainties across the sample, i.e. as 
would be introduced by assuming a uniform spectral 
shape for the stellar and dust SEDs for all galaxies be- 
tween 3.6 and 4.5 /im. 

Specifically, let us consider the emission at frequency 
v as the superposition of the emission from N different 
sources (stellar and non-stellar) such that, for P samples 
(viz. pixels) in M frequency channels, 



x = A • s 



(1) 



where x is the MxP measurement set, A is an MxN 
mixing matrix containing the coefficients of the combina- 
tion and the NxP matrix s represents the source signals, 
under the assumption that the sources S are separable 
functions of location and frequency (i.e. Sj—Sj(r) and 
A itj cx Sj(v)). 

In this way, we can extract measurements of the source 
signals Sj (each representing cither contaminants or emis- 
sion from the true old stellar disk) as some linear com- 
bination of input channels Xi (i.e. IRAC channel 1 and 
2 images at 3.6 and 4.5 fim). But rather than assume a 
particular set of separation coefficients A~- (given some 
empirical recipe), we can solve for these simultaneously, 
letting source extraction become a problem of so-called 
blind source separation (BSS). 

Several techniques have been proposed to reliably solve 

19 A recipe for scaling and subtracting the 3.6 and 4.5 iim im- 
ages to remove contaminants could be established, e.g., using a 
calibration to 8 fim data. Note, however, that simply subtract- 
ing a scaled version of the 4.5 /im image from the 3.6 fim image 
cannot account for the fact that dust and stellar SEDs have dif- 
ferent shapes between these two bands, and a single scaling is not 
physically realistic throughout the entire galaxy. 



problems of this type, including Principal Component 
Analysis (PCA) and its increasingly popular affiliate 
Independent Component Analysis (ICA). Both extract 
source measurements via transformation of the input 
data, but whereas the former imposes only that the 
sources be uncorrelated (i.e. they have zero covariance), 
the latter maximizes their statistical independence. This 
grants greater uniqueness to ICA solutions, promising 
improved descriptions of sources that are the result of 
distinct physical processes. As demonstrated in § [4j ICA 
decomposes the images at 3.6 and 4.5 /jm each into two 
other images, tracing either the old stars or the contam- 
inant sources of emission, by maximizing the distinction 
in the [3.6]-[4.5] color of these components. 
Here we employ the FastICA realization of noise-free 



ICA develo ped by |Hyvarinen| ([1999) and Hyvarinen & 
Oja| ( |2000j ) . 20 This powerfully efficient method uses a 



negentropy measure of non- "Gaussianity" to define the 
statistical independence of the sources. By the central 
limit theorem, a mixture of independent variables is more 
Gaussian than either of the two alone, and so the require- 
ment of variable-independence can be replaced with the 
equivalent requirement for less "Gaussianity" in these 
variables. 

Thorough descriptions of the method can be found, 
e.g., in the analysis of mu lti-channel COBE and WMAP 
obser vations of the CMB (| Maino et al.|2003| |Maino et al. 
2007 1 . To summarize briefly here, the separation of non- 



gaussian sources is achieved through linear combinations 
of the input data at different frequencies (e.g. y=W • x 
with y~s and W^A -1 ) corresponding to maxima of the 
negentropy. The negentropy, itself, is approximated by 
a non-linear re-mapping of the input data according to 
one of several functions g of the data's iteratively up- 
dated principal components u, the choice of which con- 
trols the form of the identified source distributions. In 
this work, we employ the fun ction g(u) cx exp(u 2 ) (see 
e.g. Hyvarinen & Oja 2000), which is thought to be 
well-suited to cases where the independent components 
are highly super-Gaussian (IHyvarinen 1999 as opposed 



to the kurtosis function, which is more sen sitive to out 
liers in the source distributions) and which |Maino et al. 
( 2007 1 find achieves the best recovery of CMB frequency 
scaling. 

By the nature of the method, sources determined with 
ICA can take on negative values since they are assumed 
to have statistically independent distributions with zero 
mean. Nevertheless, we favor this method to a third 
technique, Non-negative Matrix Factorization (NMF), 
which, despite promising physically realistic solutions 
(given that sources should be positive), is more prone 
to non-unique solutions whereby true source characteri- 
zation is less immediately obvious. 

Althoug h sparsene ss constraints on NMF (e.g. 



NMFSC; Hoycr 



2004) have been shown to improve the 



uniqueness of solutions in some applications, neither the 
sparseness of the sources nor their 'mixing' coefficients 
can be specified individually. For our purposes, then, we 
find that any such separation of sources is still less sat- 
isfactory than that which can be achieved with ICA; we 

20 We use an implementation of FastICA from within 
the IT++ library, pu blicly available for download at 
|http: / /itpp.sourceforge.net| 
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expect at least one source (i.e. emission from old stars) 
to appear both throughout the disk (in each pixel) and 
in the two wavebands. Any source negativity imposed by 
ICA is furthermore largely relegated to noise, which we 
find predominantly hosted by one source map (see eq. [2] 
below) . 

2.2. Application to S^G images 

Given two input samples (here the two images at 3.6 
and 4.5 fim) the method supplies at most two indepen- 
dent sources (i.e. M=N for A invertible), imposing a 
strict limit on the number of contaminants we can de- 
tect, in addition to emission from the old stellar disk. As 
discussed in § |5.3| however, we find that two input images 
may be sufficient for the binary distinction between old 
stars and contaminants, given that the latter tend to be 
detected together as a single second source. When this 
second source is removed from the first source, however, 
artifacts are introduced, as discussed in 



6.1 



With ICA measurements of the two sources, then, 
along with the 2x2 separation matrix W ~ A -1 , the 
components j at a given pixel in each image i can be 
reconstructed as = W j j 1 Sj (without implicit sum 
over j). (Note that for any BSS technique there is a 
scaling ambiguity between the mixing and sources, Ay 
and Sj, but the contribution of each source to the to- 
tal flux can be written as Ay-Sj without loss of informa- 
tion.) By this definition, the flux of source j in image i=l 
at 3.6/itm is related to its flux in image i=2 at 4.5 pan 
as W£ 6 j/W^ j. So, taking source si as the old stars, 
for example, this relation defines an effective corrected 
global [3.6]-[4.5] color for the si component according to 
the ratio A3. g.i/Ai. 51. 

Multiple realizations of the sources and their colors 
can be obtained by seeding the ICA calculation with a 
non-random initial guess for the mixing coefficients with 
controlled variation. Our seeds conservatively suppose 
that the flux at 3.6 and 4.5 fim is split between the stars 
and the dust such that these components have colors 
-0.2<[3.6]-[4.5] stars <0 and 0<[3.6]-[4.5] <Joat <1.5, respec- 
tively. The choice of seed is arbitrary in the sense that the 
colors of the calculated components quickly move away 
from the initial values and approach the final converged 
values, but this step allows us to define a measure of 
the uncertainty in the final solution: we take the aver- 
age and rms variation of the colors returned with our 
fixed number of seeds and then adopt the average to re- 
construct the quantity AySj for each source at the two 
wavelengths. 

2.3. Comments on ICA with S^G Data Products 

To best utilize ICA for source identification in this par- 
ticular application requires that the emission intrinsic to 
the galaxy can be defined in the input images. This pre- 
vents the ICA-identified distributions from being skewed, 
or biased, representations of the true distributions of 
emission sources. Measurements of the shape and lo- 
cation of the effective edge of the old stellar disk, such 
as measured by S 4 G pipeline III (by S 4 G convention a t 
the 26.5 magnitude/arcsec 2 isophote; Sheth et al.||2010| ), 
are therefore a crucial input for this technique. The pa- 
rameters of the edge of the disk listed in Table [T] were 
chosen by eye to enclose the disk emission according to 



the S 4 G definition. Variation by ^20% in these values 
has negligible effects on the results, whereas significant 
differences arise when, e.g. the radius is defined to be 
half as large. 

Masks of bad pixels (e.g. S 4 G masks prepared with 
SExtractor) also provide a convenient way to avoid fea- 
tures in the input images, including field stars and back- 
ground galaxies, that are clearly not inherent to the 
galaxy. However, we do not use specially designed masks 
in the ICA analysis that follows. In testing we find that 
this choice has very little effect on the results of the 
method for this sample of galaxies, despite the risk of 
bias (i.e. accurate source extraction with ICA requires 
that the intrinsic shape of the underlying source distri- 
bution be correctly identified); for the six galaxies here, 
the identified distributions vary little in cases with and 
without field stars masked. Such masking may be neces- 
sary in other cases with stronger foreground contamina- 
tion, however, and should add decided value at a stage 
of analysis beyond the application of ICA (e.g. GALFIT 
decompositions of the surface brightness distribution), 
for accurate estimation of the total flux and mass of the 
old stellar population. Other features with the potential 
to bias the results of the analysis, such as observational 
artifacts, are identified and excluded at a later stage, as 
described in § |2.4| 

Since source separation with ICA leverages color infor- 
mation, the small differences in the 3.6 and 4.5 /jm PSFs 
are a potential influence on accuracy of the solutions. 
Here we have used the convolution kernel calculated by 



Aniano et al. (2011) to match the 3.6 and 4.5 /im im- 
ages^ The results in this case are very similar to those 
obtained using unmatched images, though, with less than 
0.01 mag difference to the final source colors. The choice 
of kernel (i.e. substituting the Gordon et al. 2008 ker 



nel) or inp ut PSFs (using the S 4 G "super" -PSFs; Sheth 



et al.|20i0 ) also does not significantly impact the results. 



1'his is to be expected given that the PSFs for the 3.6 
and 4.5 pm bands are co mparably sized (a lthough highly 
structured in the wings ; |Reach et al] 12005] and [A mano 



et al. |2011 ). Any effects due to lingering PSF mismatch 



can also be minimized during the a posteriori exclusion 
described in § |2.4| For application to the full S4G sam- 
ple of extended sources the necessity of PSF-matching is 
therefore unclear, but care may be warranted in select 
cases (i.e. bright or saturated point sources in the field). 

Prior to source separation with ICA, the sky back- 
ground is subtracted from the 3.6 and 4.5 fim images for 
each galaxy. Here we define the sky level in either band 
as the average of pixels with signal-to-noise (S/N) <2 
beyond the edge of the disk. We note that the adopted 
sky level has very little influence on the identified sources 
(except for a uniform shift in flux level) , and does not af- 
fect the scaling between them. Colors listed in Table [2] 
are th erefore only subject to photometric errors (defined 
as by Salo 2010a plus ~20% for the sky error), which 
typically range from ~ 0.03 mag at (13. 6=18 mag arcsec -2 
to 0.2 mag at Li3.Q=21 mag arcsec -2 at 3.6 /im in a given 
pixel. 

According to eq. (fTl), the flux errors in both sets of 
output source images (one pair at each wavelength) can 
be expressed as a combination of the errors in the input 
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TABLE 1 
Sample 



Galaxy 


RA 


DEC 


Classification 


Lir/ L opt a 


8/24^m" 


D (Mpc) 


Redgei^ ) 


^edge 


PA ec / 9e ( ) 


NGC 1566 


04 ,l 20 m 00.4 3 


-54°56'16" 


(Rl)SAB(s)b 


0.49 


0.69 


17 


4.5 


0.996 


4 


NGC 2976 


09 h 47 m 15.4 s 


+67°54'59" 


SAB(s:)d 


0.49 


1.09 


3.6 


3.75 


0.599 


54 


NGC 3031 


09 h 55 m 33.2 s 


+69°03'55" 


SA(s)ab c 


0.08 


0.66 


3.6 


11.12 


0.55 


70 


NGC 3184 


10' 1 18 m 16.8 s 


+41°25'27" 


SA(rs)bc 


0.99 


0.90 


8.6 


3.25 


0.99 


2 


NGC 4321 


12 ,, 22 m 54.8 s 


+15°49T9" 


SAB(rs,nr)bc 


0.73 


0.81 


18 


5.62 


0.83 


68 


NGC 5194 


13' i 29 m 55.7 s 


+47°13'53" 


SAB(rs,nr)bc 


0.6 


0.82 


8.4 


5.62 


0.65 


90 



Note. — The parameters defining the edge of the disk (last three column s) are chosen by eye to encompass the 
bulk of the dis k emission at 3.6 ym. Classifi cations are from Buta et al. ( |2010 ), unless otherwise noted, and distances 
are taken from Munoz-Mateos et al. (2009a I. 



"Taken from |Kennicutt et al.| \2003\ 

'Calculated using the integrated fluxes measured by |Murioz- 

|Mateos et aT] ||2009b[ ) 

"tiAJd Classification | |de Vaucouleurs et al.|1991^ 

images, <7i and oi. At 3.6 fim (band 1), for example, 



1 



o-ii = 
and 



012 = 



1 - 10 



-0.4(d-c 3 ) frl+' 1 °- U - 4C ^) 



-^Q-0.4(cl-c2) 
1 _ 10-0.4(cl-c2) 



((7! + zHn°- 4c V 2 ) (2) 



in the first and second source images, respectively. Here, 
z is the ratio of the zeropoints in each measurement band, 
and ci and c 2 are the colors for each source (e.g. as listed 
in Table [2]). The noise in the final source images at the 
second wavelength (band 2), 021 and er 2 2, are equivalent 
to an an d C12 with only a scaling determined by the 
color of the source. In the following, photometric errors 
a m refer to these errors on the fluxes measured in a given 
source image m (dropping the waveband index). 

2.4. A posteriori exclusion 

Both source maps exhibit typically noisy fluctuations 
beyond the bulk of the emission, where the signal-to- 
noise in the original maps is low. Only rarely where the 
S/N is otherwise high will negative values persist, pri- 
marily in map s2. This 'ghosting' is symptomatic of the 
incompatibility of the true sources with the set of only 
two that we calculate here with ICA (as identified also 
by Maino et al. 2003). This could be intrinsic to sources 
in the disk (e.g. stellar distributions with distinct old 
and young components), but here it occurs mostly as 
the result of observational artifacts, i.e. in the presence 
of a saturated nucleus or due to PSF mismatch around 
bright central sources. In such cases we recalculate solu- 
tions with any such pixels masked in the input measure- 
ment sets. For saturated nuclei a circular mask is cho- 
sen by eye around the central region, while in all other 
cases masking occurs where negative- valued pixels exceed 
|10cr m | (i.e. where \S/N\>10) in either source image at 
3.6 \xm. If the number of pixels in the two recalculated 
source maps is in total reduced, these solutions replace 
the original set. Flux information from the masked re- 
gion is then passed untouched into either one of the two 
component maps, depending on which of the first source 
solutions defines the a posteriori mask. In cases where 
negative values initially appear in s2, for example, the 
new map si regains the original flux in the masked re- 



gion, while pixels in the new s2 are there set to zero. 
Note that, with a posteriori exclusion, the final maps 
reflect a combination of the uniform ICA source and a 
subset of pixels with flux and color information from the 
original uncorrected image. 

We emphasize that these instances are limited (here it 
occurs in the central regions of NGC 5194 and M81 and 
at the location of the saturated nucleus in NGC 1566, 
covering less than 1% of the disk area in each case). 
Negative-valued pixels are otherwise mostly all within 
the l<r m uncertainty for each pixel defined by Poisson 
statistics, and so can be considered equivalent to a non- 
measurement. 

In an additional post-processing step, these negative 
pixels are excluded from final ICA maps, together with 
all pixels below a 5o~ m flux threshold, where o~ m for each 
source image is calculated from the noise in the original 
3.6 and 4.5 fim images according to eq. [2j Pixels in cither 
map with fluxes below the cut-off are then replaced with 
zero (or blank), while pixels at those positions in the 
companion map are set to the value measured in the 
original image (i.e. flux in an image below the threshold 
is added back into the companion map). 

This thresholding produces final source images that 
more realistically represent the expected spatial distri- 
bution of contaminants; in the absence of real flux from 
a contaminant, ICA s2 'contaminant' images are other- 
wise filled with noise over a large area inside the disk 
edge. The final companion old stellar images in this way 
also better portray 'contaminant-free' regions. With the 
exclusion of pixels in the ICA contaminant map below 
5<72 at 3.6 fim, for instance, pixels in the old stellar map 
at the positions of these blanked pixels are returned to 
their values in the original, uncorrected maps. (In all 
cases, the 5(72 flux threshold is by far more exclusive in 
the contaminant map than the 5<7i threshold is in the 
map of the old stellar light, so that flux is more often 
reintroduced to the old stellar map, over a larger frac- 
tion of the disk.) 

This a posteriori adjustment to the source fluxes also 
varies the ratio of the sources between 3.6 and 4.5 /im 
from the global values designated by ICA. For instance, 
the [3.6]-[4.5] colors in the contaminant free regions of the 
old stellar maps take on the values originally observed in 
the uncorrected images. According to our adopted 5o2 



G 



flux threshold, this change to the old stellar color is on 
the order of the la photometric error in the two bands, 
given the typically ^20% difference in errors between the 
3.6 and 4.5 \im bands. (From equation [2j the change in 
the old stellar map color with the reintroduction of flux 
from the contaminant map below the ha-i threshold is at 
most equal to the difference between the 5er errors in the 
two bands, i.e. ha^^-ha^.^, as defined in § 2.3| Higher 
thresholds can be used at the expense of the [3.6]-[4.5] 
color.) In what follows we take the ICA [3.6]-[4.5] color 
as representative of the true old stellar disk, and assert 
that this ICA color is realistic to within la throughout 
the entire final map of the primary, old stellar source. 
Note that, without ICA, the [3.6]-[4.5] color can be diffi- 
cult to reliably measure in the outer portions of the disk, 
given the low S/N there. However, with ICA we gain 
a measure of the [3.6]-[4.5] color that is good to la, ac- 
cording to our thresholding, given the use of information 
from throughout the rest of the disk. 

3. THE S 4 G DATASET 

The analysis presented in this work is based on 3.6 
fim and 4.5 /j,m S 4 G archival images which have been 
processed through the S 4 G reducti on pipeline (first col - 



(20081. 



gale 



umn in Figure lb, as described by Sheth et al. (2010) 



The subset of galaxies considered here, JNCC 1566, JNGO 
2976, NGC 3031 (M81VNGC 3184, NGC 4321 and NGC 
5194 (M51) (see Table HI), were chosen for their range of 
morphological properties, du st contents, and star for ma- 

and 

see, e.g '' ' """ 



tion histo ries from SINGS (Kennicutt et al 

Draine et aL"1|2007[ | Calzett i||20071 |Bendo < t 



aT1[2008l [Mufioz-Mateos et al.| |2009b, Moustakas et al. 
2010|). In addition to a fairly wide range in FIR-to 



optical flux ratio (Table [IJ-which is correlated with dust 
optical depth, dus t temperature and inclination (Ken- 
nicutt et al. 2003[ |-these galaxies span a range in nu- 
ciear activity Seyfert (NGC 1566, M51), LINER (NGC 
3031, NGC 4321), a nd starburst (NGC 2976, NGC 3184) 
( Smith et al. 2007 ) . They also show a diversity in PAH 
emission relative to hot dust, as measured by the ra- 
tio 8/24 am (T able [It, t y pical for normal star-forming 
galaxies ( |Dale et al.|p005[ ). Finally, typical of SFHs for 
Hubble types Sa-Sd this sample of six galaxies covers 
2.8<FUV-[3.6]<4.5, with redder colors i ndicating lower 
current-to-pa st star formation activity (|Munoz-Mateos 



et al. 2009b). Although few in number, this subset is 
expected to be diverse enough to demonstrate the pro- 
jected effectiveness of the ICA technique over a range of 
Hubble Types. 

We compile a set of comparison data for these galax- 
ies from SI NGS (namely IRAC 8 a m and MIPS 24/xm 
images; see |Kennicutt et al. l |2003[) and 2MASS (here, 
primarily H-band images; see |Jarrett et al.||2000 ). Addi- 
tional processing of these images will be described when 
referred to in upcoming sections. 

Unless otherwise specified, all images are presented in 
units of MJy/sr. For conversion to mag/arcsec 2 relative 
to Vega, we use the zeropoint fluxes 280.9, 179.7 and 



64.13 Jy at 3.6, 4.5 and 8 pm, respectively (Reach et al 



2005 1, and the 2MASS H-ba nd magnitude zerop oint pro 
vided in the image header ( Jarrett et al.|[2000[ ). Surface 
brightnesses in units of Lq fpc 2 are also derived using the 
value of the absolute magnitude of the sun in the IRAC 
3.6 pan band, M| 6 = 3.24 mag, calculated by Oh et al. 



4. MAPS OF THE SOURCES IDENTIFIED WITH 

ICA 

Because the distributions in the source maps at 3.6 
/im and at 4.5 /im are identical to within a difference in 
scaling (denoted by the color), we present only solutions 
at 3.6 /im here. These are shown for each galaxy in the 
third and fourth columns of Figure [l] hereafter sl-ICA 
(the old stellar component) and s2-ICA (the contami- 
nants). The original 3.6 pm image is displayed together 
with the original map of [3. 6] -[4. 5] colors in first and sec- 
ond columns, respectively, and the 8 fj,m SINGS image 
is shown in the last column. 

In general, map si is a clean, smoothed version of the 
uncorrected image, while map s2 contains most, if not all, 
of the bright and knotty features tracing, e.g., the spiral 
morphology. The separation has imperfections, though, 
such as an over-subtraction of the contaminants in iso- 
lated zones that, in some cases, leaves a 'hole' in the old 
stellar light. A discussion of image imperfections will be 
postponed until § |6j 

Comparison with the [3.6]-[4.5] color maps in Figure 
[T] demonstrates that this separation of the input distri- 
bution into two distinct emission sources is effectively a 
separation by intrinsic color; regions with [3.6]-[4.5]>0 
appear exclusively in the s2-ICA map, whereas map sl- 
ICA is composed of sources with [3.6]-[4.5]<0. The final 
colors of the two sources listed in Tableware in fact con- 
sistent with either 'red' contaminant (dust) emission or 
'blue ' stellar light dominated by K or M ty pe giant stars 



(e.g. Pahre et afl[2004a Pahre et al||2004b). 



would expect for two physically distinct emission pro- 
cesses. But, even for the small sample here, each can 
be related to other observable disk properties. A larger 
sample will be essential for investigating dependencies on 
dust properties, metallicity, and star formation history, 
but we note that stellar disks assigned the bluest [3.6]- 
[4.5] colors (marked by more negative [3.6]- [4.5] colors in 
Table [2]) are those with the reddest optical colors (e.g. 
g — i MJ.85, using the asymp totic magnitudes measured 
by |Munoz-Mateos et al.|2009b ) and earliest Hubble type 
in the sample. This confirms that the ICA-measured 
[3.6]- [4.5] colors of stellar dis ks reproduce the sequence 
noted by Pahre et al| ( |2004a[ ), namely that stellar disks 
with blue [3.6J-[4.5J colors are older than their redder 
counterparts (in light of the modest age difference be- 
tween K and M giants dominating the observed light). 

In addition, the reddening of contaminants over the 
range in [3. 6] -[4. 5] color observed here coincides with 
trends in the measured proper ties of the emitting d ust in 
these galaxies, as compiled by Draine et al. (2007) using 
information from throughout the dust SED. For exam- 
ple, galaxies with the reddest dust (i.e. largest positive 
[3.6]-[4.5] in Table [2]), which have lower local contribu- 



5.2), 



tions from contaminant emission (as measured in 
tend to have smaller dust masses and a lower fraction 
of dust in the form of PAHs (measured by qpAH', see 



Old K and M giants exhibit blue colors in the IRAC bands 
due to photospheric CO absorption in the 4.5 fim band. M stars 
have colors [3.6]-[4.5]~-0.15, while younger K stars with weaker 
CO absorption app ear redder ([3.6l-[4.5l~ -0.06; e.g. as noted by 
|WilIner et al.|2004| and |Pahre et al|2004b) . 
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Fig. 1. — From left to right: S 4 G pipeline-processed IRAC channel 1 (3.6 /im) image ( jSheth et al.|2010| ); [3.6]-[4.5] color map; 
separation of sources ac hieved with ICA at 3 .6 pm, map sl-ICA (old stars) and map s2-l(JA (c ontam inants); SINGS IRAC 
channel 4 (8 /im) image ( Kennicutt et al.|2003 1 shown with the stellar continuum subtracted as in § 5,1.1| All images are shown 
on a square root intensity scale (range given below), except for the [3.6]-[4.5] color (shown on a linear scale from -0.3 to 0.3 
mag, with [3.6]-[4.5]=-0.3 in black, [3.6]-[4.5]=-0.15 in blue, [3.6]-[4.5]=0 in green, [3.6]-[4.5]=0.15 in orange and [3.6]-[4.5]=0.3 
in white). The upper end of the intensity range I max for images shown in columns 3, 4 and 5 (sl-ICA, s2-ICA, and the stellar 
continuum subtracted 8[im emission, respectively) is set by the value J™™ adopted for the original 3.6 /im image shown in 
column 1 (where 1^=5 MJy/sr in NGC 2976 and NGC 3184, I^ x =12 MJy/sr in NGC 4321 and NGC 1566, I^ x = 15 MJy/sr 
in NGC 3031 and Ig l g x =20 MJy/sr in NGC 5194). For all galaxies: IT" X =IT$ X , T& ax =Q.M¥$* and I^^lM^t , except 
for NGC 3031 where I]^l 6 =Q.2I^i x =Q.U^ ax . (In NGC 3031, where L IR /L ovt in Table JTlis lowest, contaminant fluxes are 
fainter, and the implied 3.6/8 /j,m PAH ratio lower, than in all other galaxies; see Figure [2]) Field dimensions (top to bottom): 
10.6' x 8.75' (NGC 4321); 12.26'xl3.81' (NGC 5194); 9.51'x8.96' (NGC 3184); 7.54'x6T6' (NGC 2976); 10.52'x8.82' (NGC 
1566); 22.52'x22.52' (NGC 3031). 
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TABLE 2 
ICA [3.6]-[4.5] COLORS 



Galaxy 


si 




s2 


S S 2 


NGC 1566 


-0.132 


0.021 


0.669 


0.242 


NGC 2976 


-0.143 


0.022 


0.44 


0.029 


NGC 3031 


-0.102 


0.015 


0.239 


0.194 


NGC 3184 


-0.127 


0.051 


0.764 


0.311 


NGC 4321 


-0.071 


0.019 


0.242 


0.062 


NGC 5194 


-0.094 


0.029 


0.542 


0.231 



Note. — Global [3.6]- [4.5] colors for the two ICA sep- 
arated sources, sl-ICA (old stars) and s2-ICA (contami- 
nants). Variation in the color of each component over a 
fixed range of initial seeds (see text) defines the uncer- 
tainty measure S. The larger uncertainties for source s2 
reflect the larger intrinsic spread in contaminant sources 
collected together in the single contaminant map, as in- 
ventoried in § [5j For reference, all colors wo uld be lower 
by -0.0 12 using the zeropoints adopted by |Pahre et al| 
(12004a I. 



Draine et al.||2007 Tables 4 and 5). A larger sample will 
be necessary to link these trends. In what follows we 
first undertake a detailed study of the two ICA sources 
at 3.6 fim using the spatial and intensity information in 
the ICA source maps. 

5. THE CONTAMINANTS 

In this section we explore how well the expected sources 
of contamination in the 3.6 /zm images are detected with 
ICA. These contaminants are defined in what follows as 
all sources of emission apart from the oldest stars: PAH, 
hot dust, and intermediate-age red supergiant (RSG) and 
asymptotic giant branch (AGB) stars. With the removal 
of these contaminants from S 4 G images, the remaining 
old stellar light should be readily convertable into mass 
maps. 

5.1. Emission from Dust 

The apparent similarity between maps s2 and the 8 pm, 
images for all galaxies is a first strong indication that the 
emission in the ICA-calculated map can be associated 
with non-stellar emission. A more strict identification 
of, e.g., the expected 3.3 \xm PAH contaminant there 
is complicated, however, given that the scaling between 
this feature and the standard 8 /jm PAH pr oxy depends 
on the size and ionization state of the PAH ( |Bakes et al. 
2001), and is thus fairly uncertain. 

Emission from hot dust as traced, e.g., by 24/im im- 
ages, which is present to varying degrees in these galax- 
ies (see the 8/24 /im flux ratio in Table [T]), also con- 
tributes at 3.6 [im and 4.5 /im. Hot dust emission is well- 
correlated with 8 /im PAH emission only on the largest 
scales. Below a few hundreds of parsecs (on the scale 
of HII regions) Bendo et al.| ( |2008[ ) find an excess of hot 
dust emission relative to the PAH. 

These two factors are of critical importance for inter- 
preting the s2 contaminant map. Emission from both hot 
dust and PAH (the 3.3 /im fea ture in the 3.6 /xm band 



and continuum at 4.5 /im; i.e. Flagey et al. 20061 con 
taminate the 3.6 and 4.5 /im images whereas, m perform- 
ing ICA with only 2 bands, we assume only a single non- 
stellar source. This effectively supposes that the PAH 



emission scales perfectly with the hot dust, i.e. given a 
global scaling between the 3.3 /xm PAH emission (and 
continuum) and the broad component at 8 /im, together 
with an equivalently uniform scaling between the 8 \xm 
emission and the hot dust as imaged at 24/im. 

In this case, only the dust spectral index given by 
the ratio of the s2 contaminant map at 3.6 and at 4.5 
\xm (i S 2,3.6/is2,4.5) would deviate from what we expect 
based on thermal dust emission alone. We find this to 
be the case in all our maps. The colors for source s2 
listed in Tableware bluer than the [3.6]-[4.5]~0.97 color 
characteristic of power-law emission on the Wie n side of 
th e dust SED, with spectral index rihotdust=-^ (Blain et 
al. 2003). The PAH-to-hot dust ratio is evidently close 
enough to globally constant that it can be modeled as 
such, even if there may be real local deviations (i.e. in 
HII regions). Letting ICA choose the optimal balance, 
these deviations remain recorded within our uncertainty 
measure 5 S 2, which serves to quantify the risk we in- 
herit in the limit of detecting only one contaminant where 
mor e than one source of emission may exist (see e.g. § 

This simple scenario also holds well enough that we 
find a good correlation between the contaminant flux at 
3.6 /im and the 8 /im emission, as demonstrated below. 
Even when galaxies exhibit extreme spatial variations 
in the PAH-to-hot dust ratio, and the 'perfect scaling' 
scenario becomes less appropriate, the effect on both 
the ICA-identified contaminants and old stars is small, 
within the (lcr m ) photometric errors. 

5.1.1. Comparison with 8 /im PAH emission 

Below we explore the correlation between the emission 
in map s2 and emission from dust as imaged at 8 /im, 
which is dominated by the broad PAH feature. While 
we do not expect a simple, one-to-one relation (given the 
factors discussed above, as well as the inevitable uncer- 
tainty in the subtraction of the stellar component from 
the 8 /im map), the trends are useful for distinguishing 
between hot dust, PAH emission and 'red kn ots' dom- 
inated by intermediate-age stars (e.g. § 5.2). In ad- 
dition, this exercise may motivate future initiatives to 
study PAH and hot dust emission within the s2 contam- 
inant maps; these maps of the dust have higher angular 
resolution (i.e. 1.2" in 3.6 and 4.5 /im images) than maps 
at longer wavelengths. 

Finally, comparison to the 8/xm emission provides as 
accurate quantification of the quality of our contaminant 
correction as is currently possible, given that corrections 
for the contaminants in the IRAC band s are typically 
based on such images, w here available (Kendall et al 



2008| |Helou et al.||2004[ ). Because we prefer to retain 
the high native resolution at 3.6 fim, the pixel-by-pixel 
comparisons to the 8 fxm emission are not performed on 
PSF-matched images. We have confirmed that the gen- 
eral conclusions drawn here are unchanged when the 3.6 
and 8 /im images are at the same resolution (the PSFs 
at 3.6 and 8 /xm are similarly sized, differing by a factor 
of ~1.2 in FHWM), suggesting that it is not necessary 
to degrade the 3.6 //m resolution to make identifications 
similar to those described in this section. Later in § 5.2 
though, we comment on circumstances beyond our im- 
mediate concern that require a careful treatment of PSF 
differences. 
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In Figure [2] we plot for each galaxy the flux measured 
in s2 against the non-stellar flux at 8 /xm where, for all 
galaxies, the non-stellar flux at 8 fim is calculated as 
Fs, ns = F$ — 0.232 ^1,3.6 based on th e stellar continuum 
correction given by Helou et al. ( 2004[ ), and with the con- 
vention F=f v (also used for all remaining fluxes). Note 
that in this definition we represent the stellar component 
with our map si at 3.6 \xm. 

Most galaxies ( with t he exception of NGC 3031 and 
NGC 1566; see § 5.1.2) show a clear relation between 
the flux in s2 and the PAH emission at 8 /im over al- 
most the full range in Fs. ns - Notably, the slope of the 
implied linear relation for each falls within the range of 
expected values for the scaling a between the PAH emis- 
sion at 3.3/im and at 8 jitm, 0.03<a<0.09 , as measured 
in the Milky Way (see Flagey et al. 2006 and references 
therein). At the very least, this suggests that the ICA- 
idcntified component is a reasonable map of the PAH 
emission at 3.6 fim, by virtue of the realistic fluxes as- 
signed to s2 (i.e. the ICA 'measured' fluxes are consis- 
tent with a dust component). There is obvious scatter, 
though, and the overall relation in some cases features a 
shifting slope with increasing or decreasing F S ns . 

The scatter in the correlation might be partly at- 
tributed to differences in the size and ionization state 
of the PAHs throughout the disk. But in addition to 
physical variations in the PAH, uncertainty in the stel- 
lar continuum subtraction (i.e. given error in si, or in 
the scaling of the old stellar light from 3.6 to 8 /im), a 
non-dust origin for some of the 'contaminant' flux in s2, 
or thermal dust emission could plausibly contribute to 
the overall scatter in the panels of Figure [2] The first 
case does not seem to be a major source for this small 
sample; we find no significant change in the scatter or 
overall appearance of these plots with -F51.3.6 replaced by 
the original, uncorrected fluxes in the definition of F% ins 
above. In the second case, the signature of a non-dust 
source of contamination is more obviously manifest to the 
left in the plots, where F s2j 3.6 is in exce ss above F 8 ^ ns , as 
considered further in the section § |5.2| 

The more likely contributor to the plots in Figure 2] 
where i 7 s,ns^2 MJy/sr (i.e. any trend for a variable 
slope, but also some scatter), and especially to the over- 
all appearance in the case of NGC 3031, would seem to 
be the non-PAH (hot dust) component in both the 8 fim 
and 3.6 fim bands, as demonstrated in the following sec- 
tion. 

Another possible source of correlation between -^2.3.6 
and F% >ns could arise with our definition of Fs. ns , al- 
though we do not find strong evidence for such a pos- 
sibility here. Specifically, using F s i 3 e rather than F 3 6 
to subtract the stellar continuum contribution in Fg in- 
troduces in F^^ns an explicit dependence on -^2,3.6 (i- e - 
F 8 , ns =F 8 -0.232> 3 . 6 +0.232.F S 2,3.6), potentially leading to 
a forced correlation at fixed observed F s and F 3 . 6 . How- 
ever, the latter two observables are far from fixed, and 
the variance in the quantity ^-0.232^3.6 is at least an 
order of magnitude larger than the variance in ^2,3.6 
in all cases. (In one example, NGC 4321, the variance 
in ^-0.232^3. 6 is 13.6 MJy/sr, while the variance in 
-Ps2,3.6 is 0.6 MJy/sr.) Meanwhile, the variance in the 
constructed Fg^ ns is comparable to the variance in F$- 
0.232F 3 . 6 (e.g. in NGC 4321 the variance in F 8 ,„ s is 13.4 
MJy/sr), which strongly indicates that F$ tns is driven 



primarily by Fg and F3.6, not -£^2,3.6 ■ Likewise, we see 
very little change in the appearance of these plots with 
Fsi,3.6 replaced by F3.6 in the definition of Fs,ns, as men- 
tioned earlier. Altogether, this supports the conclusion 
that Fg^ns depends only negligibly on i^2,3.6i or that the 
correlations between -F S 2,3.6 and the imperfectly indepen- 
dent variable Fg^ ns are not the result of the (minor) de- 
pendence of Fg iTls on -F S 2,3.6- Instead, these correlations 
arguably arise from the physical characteristics of the 
emission contained in map s2. 

5.1.2. Contribution from hot dust 

Comparison with the non-stellar emission at 8/im 
demonstrates not only the detection of PAH emission in 
the contaminant maps at 3.6 fim but also a second, hot 
dust contaminant. Consider, for instance, the grouping 
of pixels at low 8 fim flux (F s <3 MJy sr^ 1 ) in either 
NGC 1566 or NGC 3031 following a steepened trend. In 
both cases, these pixels can be identified with regions 
in the disk where the ratio Fs/Fm drops, and hot dust 
dominates the emission from PAHs. In NGC 1566, for 
example, these pixels cluster in small isolated HII regions 
that string along the spiral arms, down in to the center 
of the disk (see Figure |3|). In NGC 3031, these pixels are 
located in the central 75", where -Fg/i*^ is low in the 
surrounds o f NGC 3031's AGN (see Figure[3) e.g. |Bendo 
et al.|2008 1 . In the immediate vicinity of the unobscured 
nucleus in this particular case, the mid-IR emission can 
be associated with the signature "big red bump" emission 



from a truncated accretion disk ( |Ho et al.|1996 Ho|2008 



Quataert et al.|1999[ ), but emisson from an extended hot 
dust component as far out as 75" (possibly heated by the 
AGN) presumably also contributes to the drop in Fg /F24 
here. 

In these two cases we expect that the steepened cor- 
relations in Figure [2] are dictated less by the PAH than 
by the other main source of mid-IR emission, given the 
observed low level non-stellar F 8 /F 2 4. In this case, the 
ratio F S 2 j 3.6/Fg ns at each pixel should roughly measure 
the non-stellar spectral index n (but see the upcoming 
discussion). Using the measured ratio F,,2. 3.6/F 8 . ras we 
find a spectral index consistent with n=2 (Blain et al 



2003| as expected for power-law emission on the Wien 
side of the dust spectrum) with which we recover the 
observed flux at 24/xm within 15%, on average, in both 
NGC 1566 and NGC 3031. 

A large spread in the flux of hot dust over the disk 
will result in a 'continuum' of shifted linear relations be- 
tween F S 2 and F 8ins for which a may also vary, possibly 
contributing to the overall scatter in these plots. But 
because the fall-off rate of hot dust emission between 8 
and 3.6 fim is smaller than ratio of the PAH features in 
these bands, the signature of the hot dust will in gen- 
eral appear systematically to the left of the PAH in the 
F 3 6 -F 8 space plotted in Figure[2] making this component 
straightforward to identify in the contaminant map. As 
depicted in Figure [3j the range 0.1<F s2: 3.6/F 8 <0.3 can 
be used to identify pixels with emission from hot dust, 
while the criterion F s2i 3. 6^0.09 selects for PAH emission. 

The spatial locations of hot-dust dominated regions 
can be successfully identified in this way. Hot dust fluxes 
may be slightly less reliable without correction, though, 
as they are susceptible to overestimation. Since the spa- 
tially predominant PAH component often more heavily 
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Fig. 2. — Contour plots of F S 2.3.e vs. Fs, ns (calculated as described in the text) for NGC 1566 (top left), NGC 2976 (top 
right), NGC 3031 (middle left), NGC 3184 (middle right), NGC 4321 (bottom left) and NGC 5194 (bottom right), with number 
of image pixels increasing from the outermost (orange) contour to the innermost (purple) contour. Over-layed dashed lines 
represent the expected range in F S 2,3.e~aFs ns for PAH emission, with «=0.09 or 0.03. A third solid line shows the relation 
given by the best-fit slope a listed in Table [3] 



controls the contaminant color than the hot dust, then 
the ICA [3.6]- [4.5] color is slightly too blue in areas with 
the highest contribution from this source (i.e. in HII re- 
gions). 22 This leads to an overestimation in the flux from 
the hot dust there, which would nominally appear more 
red in the [3.6]-[4.5] color. But as long as the flux in hot 

22 Whereas hot dust should have [3.6]-[4.5]~0.9, for the PAH we 
expect [3.6]-[4.5]~0.04-0.2 according to the measured PAH contr i- 
butions in the 3.6 and 4.5 bands (see, e.g., |Flagey et ah (2006)) 
and as is typical for the [3.6]-[4.5] color maps in Figure [il 



dust F^fP is much larger than the flux from PAHs, the 
fractional overestimation in F^ 1 ® will be independent of 
•Fjjjg and proportional to the color difference between 
the ICA value and the true value, leaving a well-defined 
trend in 3.6-8 \im space identifiable with hot dust emis- 
sion. 

Conversely, as discussed more thoroughly in § |6.1| the 
presence of hot dust can contribute to the underestima- 
tion of the PAH emission throughout the rest of the disk, 
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TABLE 3 

PAH Contamination between 3.6 and 8 /im 



Galaxy 



fit 



NGC 1566 0.065±0.0001 F 8 „ s >5 MJy/sr 

NGC 2976 0.068±0.0002 F S ns >l.5 MJy/sr 

NGC 3031 0.036±0.0001 L 8 „ s >2.5 MJy/sr 

NGC 3184 0.049±0.0001 F$ „ s >2.5 MJy/sr 

NGC 4321 0.062±0.0002 F Sns >2 MJy/sr 

NGC 5194 0.031±0.0002 r<225", F s ns >3 MJy/sr 



Note. 



F, 



s2,3.6 : 



=<xF 8 



Parameters of the best fit straight line 
in the scatter plots shown in Figure 



with selection criteria for the fit regions listed on tne 
right. These fitted regions are chosen by eye to avoid 



5.2 ) in order 



emission from hot dust and red objects (see \ 
to most fairly portray the behavior of the PAH emission 
between 3.6 and 8 fim, expected to have Qf=-F S 2,3.6/-p8,ns 
in the range 0.03<a<0.09. 



but only most strongly in cases where the ICA [3.6]-[4.5] 
color nears that of the hot dust. The color of the intrinsi- 
cally bluer PAH emission would then be over-predicted, 
resulting in an underestimation of the flux from this com- 
ponent. 

Even in cases where the hot dust becomes more dom- 
inant overall here we find little evidence that the PAH 
is affected, and the picture of the dust presented by the 
contaminant maps remains physically realistic. Galaxies 
with the reddest contaminant color are not those with 
the lowest a, for instance, although we note that the co- 
incidence of a low overall a with evidence for a second 
(e.g. hot dust) contribution, for instance, would still be 
arguably consistent with a link between t he ionization 
state of the PAH (as indicated by a; see Bakes et al. 
2001 1 and the amount of 24^m emission, with dust tem- 
perature the commonality. Plus, at least in this small 



sample, the ICA contaminant flux is arguably the fairest 
balance between the PAH and the hot dust in these maps 
where the latter component is most prominent. Emission 
in NGC 3031 is divided into almost equal subsets of PAH 
and hot dust, for example, whereas the majority of emis- 
sion in NGC 4321 is consistent with PAH emission (see 
Figure [2}. 

Based on the demonstrated correlation with the 8 \xra 
emission, we argue that map s2 provides a realistic de- 
scription for the state of the dust in the studied galaxies. 
Uncertainties in the quality of the resulting correction to 
3.6 fim images should be kept in mind, however, given the 
modest drawback to detecting the dust features together 
as n oted earlier in this section (and addressed again in § 
5.31. 



5.2. RSG and AGB stars 

Next we reveal an additional source of contamination 
in ICA map s2 besides the combination of hot dust 
and the 3.3 fim PAH emission, namely clusters that 
host bright but low M/L evolved AGB and RSG stars 
contributing to the emission at this wavelength. Such 
evolved objects are typically clearly identifiable by their 



location in [3. 6] vs. [3.6]-[4.5 



(CMDs) (c.f. Jackson et al. 



color-magnitude dia grams 



'red' [3.6]-[4.5]>0 colors. 23 In the ICA-corrected maps, 
on the other hand, colors are at a constant value [3.6]- 
[4.5]=-0.1±0.03 on average (see Table [2| , comparable t o 



what is expected for an M giant ( jPaEre et al||2004a ) 
By extension, if regions dominated by evolved RSG and 
AGB stars are present in the original 3.6 fim image, they 
can be identified with ICA as part of the 'contaminant' 
distribution in the s2 image. 

Qualitative comparisons between the contaminant map 
and the 8 fim image are often alone suggestive of a dis- 
tinct secondary component. This is especially obvious for 
NGC 3031, where the ICA map s2 exhibits flux in loca- 
tions distinctly apart from the (strongly) emitting zones 
in the 8 fim image. Below we introduce a more quantita- 
tive characterization of the second component in the ICA 
contaminant distribution by combining the results of the 
previous section with spatial and color information. 

To start, pixels originating in regions with little to no 
8 fim dust emission can be easily identified in the plots 
of Figure [2] as a thin, almost vertical concentration of 
points at the far left (Fs }Tls <2 MJy/sr). These can be 
isolated using a coarse selection F s % t 3.6^0. 3Fg, avoiding 
the zone in 3.6-8/xra space populated by hot dust near 
i 7 s2,3.6/-Fs~0.2 (for power-law emission on the Wien side 
of the du st spectrum with index n—2, following |Blain et] 



al. 2003 as discussed in the previous section). The lo- 
cations of a representative set of these pixels (i.e. only 
those above F S 2,3.6~0.2 MJy/sr, for the purposes of clar- 
ity) for all galaxies are shown in black in Figure [3] To 
minimize contamination by obvious foreground stars, the 
selection includes an upper limit on the flux at the posi- 
tion of candidate pixels in the original image. 

The spatial configuration of the black pixels in these 
maps is clearly distinct from the dust, represented in light 
gray (PAH) and dark gray (hot dust), but in some cases 
similar to what we might expect for a random distribu- 
tion of field stars and background galaxies. (Recall that 
no masking of such objects has been invoked.) While 
external contaminants should in large have colors [3.6]- 
[4.5]~0, we find mostly colors 0<[3.6]-[4.5]<1.0 for the 
second distribution of contaminant pixels. 

This is clear in Figure [4] showing representative color- 
magnitude diagrams for the dust and 'red knot' compo- 
nents in s2-separated as described above-which we con- 
struct using the original color and magnitude at each 
pixel position in the uncorrected images (see the second 
column in Figure [I]). In each plot, the two distributions 
populate different, yet overlapping regions in this color- 
magnitude space: pixels with F S 2,3.6>0.3Fg extend to 
fainter magnitudes and have redder [3.6]-[4.5] colors (at 

a given magnitude) than pixels with J 7I S 2.3.6< 0-3J 7 8- 

According to typical CMP identifications (|Jackson et| 



|al.||2007 Davidgc 2009) the red colors of the former set 
are consistent with evolved, post-MS stars (and with pre- 
dictions for AGB stars, in particular, the mo st luminous 



with 0.5<[3.6]-[4.5]<l; |Groenewege"n~||2006[ ). These re 



23 Dusty, mass-losing RSGs will have colors in excess of the blue 
[3.6]-[4.5] color exhibited by other stars of K and/or M spectral 
type, depending on metallicity Z (mass loss and dust production 
increase with increasing Z). Even the youngest of these may then 
appear within the contaminant map where the effect of the dust in 
the IRAC bands exceeds the CO a bsorption at 4.5 um (which also 
increases wi th increasing Z) . (See |Bonanos et al.| | |2010[ l , |Jackson"| 



their characteristically (though not necessarily uniform) 



2007[ |Davidge|[2009[ ) given [etall {2007| l and references therein} 
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Fig. 3. — Map of pixels representative of the 'dust' and 'red knot' components in the ICA contaminant distributions for NGC 
1566 (top left), NGC 2976 (top right), NGC 3031 (middle left), NGC 3184 (middle right), NGC 4321 (bottom left) and NGC 
5194 (bottom right; only the inner r<225" is shown). Emission from PAH (light gray), hot dust (dark gray) and 'red knots' 
(black) is selected from the range F S 2,3.6/^8,ns<0.09, 0.1<F S 2,3.6/Js,ns<0.3, and i ? S 2,3.6/fs,ns>0.3, respecitvely. 
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both PAH and hot dust) and secondary 'red knot' (F S 2.3.6>0.3Fs,„ s ; in black dashed line) components in the ICA contaminant 
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(bottom left) and NGC 5194 (bottom right). (See text for selection criteria.) 
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gions could also be sites of very recently formed massive 
star clusters prevented from d ispersal that illumin ate the 
surrounding diffuse dust (e.g. Hunter et al.|2006 ) and so 
appear red in the IRAC bands. A subset are also likely 
identifiable with nearby background galaxies for which 
the obser ved mid-IR colors are closest to rest-frame col- 
ors (e.g. Fazio et al. 2004). These possibilities will be 



considered in greater detail in future work using archival 
broadband data to construct the optical-to- mid-IR SEDs 
in these regions (Meidt et al. , in prep) . Here we note that 
preliminary analysis in NGC 4321 suggests these objects 
are intrinsi cally red, with optica l colors consistent with 
AGB stars ( |Jackson et"al~1|2007| V-I -1.5). In addition, 
the [3.6]- [8] colors in these regions are indicative of emis- 
sion from circumstellar dust, as expected around A GB 
stars undergoing high mass loss ( | Jackson et al.|[2007| ). 

In any case, it is clear that this second contaminant 
population is distinct from the greater bulk of PAH and 
hot dust emitting at 3.6 /im as depicted in Figure [3j and 
that a set of relatively young stars (with ages between 30 
Myr and 1 Gyr) dominates the emission. Note, though, 
that the maps in Figure |] may strongly reflect a selec- 
tion bias against spatially coincident dust and clusters 
dominated by intermediate-age stars. At the extreme, 
the redder component in the contaminant map for NGC 
2976, for instance, is found preferentially outside the bulk 
of the dust emission. This likely arises with our iden- 
tification of red clusters in 3.6-8 fj,m space (Figure [2]), 
which requires low values of Fg ifls (i.e. little emission 
from other dust sources) to detect pixels with the high- 
est F S 2^ s .e/ Fg,^ns- A more sophisticated selection would 
be required to map this type of pixel where, e.g., F S 2 iS .6 
is high and Fs tTls is not simultaneously low. 

Despite this over-simplicity, the separation criteria de- 
veloped here are suitable for obtaining rough estimates 
of the individual contributions of the two contaminant 
components to the observed emission, as listed in Ta- 
ble |U Entries there indicate the fraction of contaminant 
emission per pixel averaged over N pixels, where N is 
the number of contaminant pixels in the case of 'local' 
measurements (columns 3 and 5), or the number of all 
pixels in the image out to the radius given in Table [T] 
in the case of 'total' (columns 2 and 4). Note that the 
derived fractions are influenced by location-dependent ef- 
fects, i.e. a reduced contribution from old stellar light at 
larger radii. 24 In addition, although averaging should 
reduce the effect of lingering cross-contamination in the 
two components, an additional bias may stem from a lack 
of PSF-matching between the 3.6 and 8 \im images; the 
8/1771 comparison directs our association of contaminant 
pixels with one component or the other. Each measured 
fraction also inherits the risk that more than one source 
contaminates any given pixel. We quantify this using the 
uncertainty measure for each contaminant color listed in 
Table [2j assuming that the color uncertainty is driven 
entirely by the 3.6 /im flux; for each fraction /, the error 
df=fS S 2 by this definition. 

As expected, we find that dust can contribute between 

24 Given the underlying radial decrease in the old stellar light, 
the pixel positions mapped in Figure [3] must be taken into consid- 
eration when comparing the 'local' values for the two components. 
(In NGC 2976, for example, the contribution of the old stellar light 
to pixels with emission from AGB/RSG stars is much lower than 
for pixels with dust emission.) 



~3 to 10% of the total emission imaged at 3.6 /im, while 
light from intermediate-age, post-MS stars accounts for 
only Locally, however, the contribution from ei- 

ther contaminant to the observed surface brightness can 
be much higher. We note that the local measurements 
in Table [4] in particular may be biased high; as com- 
mented on later in § [6j the contribution in regions with 
intrinsic colors much in excess of the ICA contaminant 
color (e.g. those dominated by AGB/RSG stars) will 
tend to be overestimated in the contaminant map. The 
dominant old stellar light keeps this fractions accurate 
to within ~10%, though. In addition, variation in the 
level of contamination from galaxy to galaxy may arise 
with the different distant-dependent physical scales be- 
ing sampled, especially for the point-like contribution 
from intermediate-age clusters. A more careful evalu- 
ation of the contaminants is beyond the scope of the 
present work. Currently, it is clear that 3.6 [im emission 
from the PAH, hot dust, or red cluster components can 
be individually studied, and with the additional bene- 
fit of PSF-matching or aperture photometry for accurate 
flux measurement. 

5.3. On the identification of multiple distributions 

In the previous two sections we clearly identified two 
main sources of contamination in maps s2, dust and 
intermediate-age, post main sequence stars. Note the 
implication, in this case, that the contaminants are de- 
tected together as a single, second source. Although this 
is not necessarily guaranteed by the method, the underly- 
ing distributions of the emission components in this par- 
ticular application offer some insight into such a scenario. 
Consider that the fraction of pixels containing emission 
from the old stellar disk is much larger than the fraction 
of pixels containing significant emission from a contami- 
nant. In addition, we might expect the contaminants to 
be more sparsely defined within a given image, as well as 
possibly from image to image, than the old stars. The 
combination of these two conditions would seem to make 
the distribution of old stars more readily definable than 
the distribution of contaminants, which then follows de 
facto. 

Only with a third measurement sample (e.g. another 
input image) is it possible to separate s2-ICA into two in- 
dependent sources (dust emission and AGB/RSG stars; 
for three sources in total), and therefore unambiguously 
identify a third unique component. However, any addi- 
tional separation would not likely grant much improve- 
ment in the corrected maps over what we derive here, 
especially since a third band will trace a slightly differ- 
ent set of (contaminant) sources than those appearing in 
the 3.6-4.5 [im window. 

Moreover, as long as the combined contaminant 
distribution-which must be more Gaussian than cither 
of the two statistically independent distributions-is itself 
still significantly different from Gaussian, then it can be 
identified with ICA, and so can its counterpart (here the 
corrected image sl-ICA representing the old stars). 25 In 
addition, as discussed below, we find that the detection 
and removal of contaminants using only two bands has 

25 Note that this explicitly assumes that our maps identified 
as s2-ICA, and not sl-ICA, reflect more than one source, as the 
preliminary comparisons above would seem to suggest. 
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TABLE 4 

Emission contributed by contaminants 



Galaxy 


s2, dust 
total 


local 


s2, stars 
total 


local 


NGC 1566 
NGC 2976 
NGC 3031 
NGC 3184 


0.1±0.024 
0.17±0.005 
0.05±0.009 
0.03±0.009 


0.34±0.082 
0.28±0.008 
0.08±0.016 
0.19±0.059 


0.01±0.002 
0.01±0.0003 
0.001±0.0002 
0.007±0.002 


0.7±0.17 
0.59±0.017 
0.43±0.083 
0.24±0.075 



NGC 4321 0.1±0.006 0.29±0.018 0.01±0.0006 0.84±0.052 
NGC 5194, r<225" 0.08±0.018 0.21±0.048 0.001±0.0002 0.38±0.088 

NOTE. — Average fraction of the 'total' (columns 2 and 4) and 'local' (columns 3 and 5) emission contributed by the two 
components in map s2 representing dust and intermediate-age RSG and/or AGB stars, with 'total' and 'local' as defined in the 
text. The stellar (dust) component here is identified with the selection i ? S 2,3.6>0.3i ; 8 ins (F S 2.3s <0.3Fs,na), similar to Figures|4] 
and [3] Errors are defined as given in the text. 



a very mi nima l effect on the corresponding old stellar 
map (see 



5.1) 



6. QUALITY ASSESSMENT OF ICA SEPARATION 

The two ICA sources, detected via combination of the 
two input measurement bands at 3.6 and 4.5 fJLm, display 
remarkable consistency with the expected characteristics 
of the main sources of emission at 3.6 /jm. Removal of 
the contaminant emission, which is well correlated with 
the non-stellar emission imaged at 8 /im, reveals an un- 
derlying old stellar disk with [3.6]- [4.5] colors that are 
consist ent with a populati on dominated by old K and M 
giants ( Pahre et al||2004"a see Table |2|. This is an espe- 
cially favorable outcome since we expect little to no dif- 
ference in the colors of stars at these wavelengths, so far 
along the Ray leigh- Jeans tail of the stellar SED. The fact 
that the colors in the final si maps are consistent with 
the colors of late-type giants suggests that each reflects 
a (nearly) uniformly old stellar population free of con- 
tamination from dust and bright intermediate-age stars. 
Final source images are not without their imperfections, 
though, as discussed below. 

6.1. Image imperfections 

Imperfections in final source images become apparent 
when the number of true sources becomes manifestly in- 
compatible with the set of only two detected with ICA 
(e.g. when the 'perfect scaling' of c ontaminants between 
3.6 and 4.5 \im discussed in § 5.1 breaks down). Large 
variations in the hot dust-to-PAH ratio, in particular, 
can reduce the effectiveness of the ICA representation of 
some part of the mapped non-stellar emission and lead 
to either the over- or under-subtraction of the contami- 
nant emission from the old stellar light. As described in 



5.1.2 a contaminant color dominated more by the hot 



dust than the PAH may result in an underestimation of 
the flux from PAH away from, e.g., HII regions with the 
highest contribution from hot dust. 

At least for this small sample any underestimation of 
fluxes appear modest, and the influence of any error 
in the contaminant color on the old stellar color even 
smaller; a change in the contaminant color introduces a 
change in the old stellar color proportional to the ratio 
of the contaminant-to-old star flux (see eq. [I]). In fact, 
adjustments to the old stellar color are on the order of, or 
within, the photometric color errors even with a boost in 



the contaminant flux by the amount implied by the range 
in the s2 color (i.e. the uncertainty measure in Tabic pi). 

Galaxies with the greatest variation in their hot dust- 
to-PAH ratio will be the most susceptible to this er- 
ror. The more frequent the HII regions, or the more 
prominent the AGN-i.e. where the hot dust becomes 
dominant-the greater the risk that the remaining PAH 
contaminant flux is underestimated. (Note, though, that 
the resulting error in the old stellar color is expected to 
remain small, given the proportionality discussed above.) 

Final source images are also susceptible to the opposite 
trend, i.e. the overestimation of the hot-dust flux when 
the PAH dictates the color of the contaminant map. This 
overestimation will be limited to small, isolated (HII) re- 
gions (at least if PAH is dominant and spatially extensive 
enough that it does indeed control the s2 color, e.g., in 
M51). Evidence for overestimation in the flux of hot dust 
is clear in Figure [2] for M51, for instance, which features 
a clear upturn at the highest values of Fg <ns , but in only 
a small number of pixels. In this galaxy and in NGC 
4321, the overestimation is large enough that t here is no 
remaining flux at the sites of HII regions (cf. |Bendo et] 
|al.||2008[ ) in map s2 (i.e. with overestimation of the hot- 
dust flux, all emission in the original maps is associated 
with contaminants). Likewise, fluxes of the RSG/ AGB 
stars with colors much in excess of the ICA contaminant 
color will tend to be overestimated in the contaminant 
map, occasionally leaving over-subtracted holes in the 
old stellar light at these locations. 

Although this is a clear loss of information in the old 
stellar maps in these regions (which would be masked, in 
any case, in the uncorrected images), the fluxes and col- 
ors elsewhere should be more realistic than given a red- 
der ICA contaminant color, since the contaminant flux is 
overall better representative of the spatially predominant 
PAH emission. 

Straightforward corrections to the old stellar map at 
the locations of these over-subtracted holes are possi- 
ble. An obvious option is to replace the holes with fluxes 
estimated from neighboring pixels, or to use models of 
the stellar surface brightness (derived with these regions 
masked) for the substitution. In both cases, the replace- 
ment should be fairer than would be possible, e.g., at 
the sites of HII regions in the uncorrected images, where 
neighboring regions may themselves reflect a non-neglible 
contribution from contaminants. 

These substitutions could then also be used to improve 
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the ICA estimates for either the non-stellar flux in HII 
regions or the contribution of intermediate-age stars to 
regions dominated by these objects, when overestimated 
in contaminant maps at the location of the old stellar 
holes. As a revision to the upper bound on the flux from 
these components in the contaminant map, an estimate 
for the flux can be established with the residuals between 
contaminant-free models and the original 3.6 /im emis- 
sion in these zones. Alternatively, it is possible to de- 
rive a correction factor for, e.g., hot dust fluxes directly 
from information in the contaminant map (i.e. the color 
and an underlying level of PAH emission estimated from 
neighboring regions), together with an assumed [3.6]-[4.5] 
color or spectral index. 

As the zero-valued holes are limited to small, isolated 
regions, we argue that this cost is worth the overall high 
quality with which the old stellar light can be mapped. In 
a forthcoming paper, we will expand the investigation of 
the ICA-identified stellar light distribution initiated be- 
low, considering, in particular, the implications for stellar 
mass maps. We find that our removal of contaminants 
enhances 3.6 fim images as tracers of the stellar mass 
surface density. 

6.2. The Old Stellar Light: comparison with NIR 

As a demonstration of the quality with which the cor- 
rected fluxes at 3.6 /im represent emission from the old 
stellar disk, here we contrast map si (and the uncor- 
rected image) for each galaxy with the H-band-observed 
distribution of light. Table [5] lists the zero-order best-fit 
to the H-[3.6] color profiles plotted in Figure [5] for the 
original and corrected maps, both of which have been 
gaussian-smoothed and sampled at the resolution of the 
2MASS H-band image. (The outer limit of the fit region 
is determined by the H-band sensitivity.) 

The bluer colors for the contaminant-corrected maps 
in all cases reflect a reduction in the amount of light 
relative to the original maps; contaminants globally in- 
troduce 0.2 mag on average across the sample. In 2D 
maps of the H-[3.6] color this reduction occurs in local- 
ized zones, i.e. in the spiral arms, and appears in the 
ID profiles in Figure [5] as a flattening of deviations from 
the mean in the strongest cases (e.g. NGC 1566 and 
NGC 5194). Notably, with emission from contaminants 
removed, maps si show greater consistency with a pre- 
diction from synthesized stellar SED fits to an old, dust- 
free elliptical for H-[3.6]~0.4, 26 although this is an im- 
perfectly matched benchmark; the current sample draws 
almost exclusively from the later Hubble types (i.e. star 
forming disks with dust emission) and so the H-band im- 
age itself may not be free of extinction and emission from 
young stars. Star formation history and metallicity will 
themselves also contribute to real differences from the 
color expected for an old elliptical galaxy. 

Metallicity may very well be contributing to the H-[3.6] 
colors here, as evidenced by radial color gradients that 
are largely positive. This signature of decreasing metal- 
licity with increasing radius has b een previously ide nti- 
fie d in these bands fo r ellipticals by Temi et al. ( 2008 1 . If, 
as Temi et al. ( 2008 1 argue in this case, enhanced, repro- 
cessed stellar emission appears in NIR and not mid- IR 

26 see |http: / / web. ipac.caltech.edu/staff/jarrett /| 
irac/calibratioii/galaxics.html 



TABLE 5 
H-[3.6] COLORS 



Galaxy Original sl-ICA 



NGC 1566 
NGC 2976 
NGC 3031 
NGC 3184 
NGC 4321 
NGC 5194 



0.81±0.004 
0.62±0.003 
0.69±0.002 
0.79±0.003 
0.77±0.003 
0.66±0.004 



0.40±0.004 
0.44±0.008 
0.54±0.001 
0.66±0.003 
0.61±0.002 
0.45±0.001 



Note. — Radially averaged H-[3.6] colors for the orig- 
inal and ICA-corrected images calculated from the pro- 
files shown in Figure [5] (Errors are statistical.) 



images, this compromises the NIR bands as abundance- 
independent measures of local stellar mass. Reddening 
associated with localized extinction (e.g. in dusty spi- 
ral arms) is also arguably evident here. Even after the 
contaminants are removed, the profiles for the corrected 
maps show lingering bumps toward higher H-[3.6]. These 
two factors suggest that map si may be a cleaner direct 
tracer of the old stellar disk than the H-band image; the 
ICA correction allows 3.6 [im imaging to better fulfill 
its potential as the optimal window on the oldest stars. 
But, to be fair, the value of ICA correction is less clear 
globally or in ID than in 2D. Non-stellar emission con- 
tributes only very modestly to the integrated light at 3.6 
y«ra and therefore amounts to a small correction in the 
global H-[3.6] color. The higher local contribution from 
contaminants tabulated in Table |4j on the other hand, 
implies as much as 0.45 mag correction to the H-[3.6] 
color in dusty regions and on average 0.8 mag change 
to the H-[3.6] color of red clusters. This clearly sug- 
gests that errors in the mass surface density of this size 
or larger (depending on the M/L) would be introduced 
with conversion to stellar mass without prior removal of 
the non-stellar emission. 

In light of this evidence, ICA-corrected images are fa- 
vorable to our end goal of mapping 2D stellar mass distri- 
butions. We note that this marks a significant advance in 
the use of 3.6 /xm images for this purpose. The presence 
of contaminants has typically involved the loss of either 
2D or structural (i.e. spiral arm) information, given the 
use of models (with an axisymmetrically-averaged spi- 



emission (i.e. 


de Blok et al. 


2008 Leroy et al. 


2008 


Kendall et al. 


2008| Kendall el 


; ai.||20ilD. 



By comparison, the independence of ICA from any 
model-inherited assumptions (including those for a par- 
ticular stellar surface brightness profile, disk shape, or 
3.3/8 fim PAH ratio) is a keen strength. Whereas cor- 
rections that rely on models are even susceptible to un- 
certainties in the assumed center position (e.g. the cor- 
rection developed by Kendall et al. 2008 that employs 
a model of the stellar continuum at 8 /xm)-which, more- 
over, need not be the same between the distributions of 
light and mass-but no such prior is required to perform 
ICA. 

Since the ICA correction can maintain full and realis- 
tic structural information at 3.6 [im (examined in greater 
detail in Paper II) it is a strong alternative to typical cor- 
rections for non-stellar contaminants in this band. Fur- 
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Fig. 5.— H-[3.6] color profiles for NGC 1566 (top left), NGC 2976 (top right), NGC 3031 (middle left), NGC 3184 (middle 
right), NGC 4321 (bottom left) and NGC 5194 (bottom right). The profiles for the original and corrected maps are shown in 
thin red and thick black solid lines, respectively, with dashed errors calculated from the range of values at each radius. The 
H-[3.6] color expected for an old, dust free population, H-[3.6]=0.4 (see text), is shown as a gray horizontal line. 



thermore, the method is simple and fast, and can be run 
largely autonomously for a sample as large as S 4 G. 

7. SUMMARY AND CONCLUSIONS 

In this paper we demonstrate the isolation of old stel- 
lar light in S 4 G 3.6 and 4.5 fxm images from contami- 
nant emission using an ICA technique designed to iden- 
tify statistically independent source distributions. ICA 
decomposes each image at 3.6 and 4.5 /im into two other 
images, tracing either the old stars or the contaminant 
sources of emission, by maximizing the distinction in the 
[3.6]-[4.5] color of these components. We find that the 
separation achieved with ICA yields a very satisfying 
view of the old stellar light that, in all cases in the small 
sample of galaxies analyzed here, has [3.6]-[4.5] colors 



consistent with the colors of K and M giants, retains a 
high degree of structural information, and fares as well 
as, if not better than, other NIR tracers as an optimal 
tracer of relative or absolute stellar mass. 

The quality of our maps of the old stellar light stems 
from the ability of ICA to detect all sources of emission 
apart from the oldest stars: hot dust, PAH and bright 
but low M/L intermediate-age stars. By comparing the 
ICA contaminant maps at 3.6 fim to the non-stellar emis- 
sion imaged at 8 \xm all three of these sources are clearly 
identifiable. In all cases we find that the ICA contami- 
nant map presents a realistic picture of the state of the 
dust and report 1) the detection of hot dust and 3.3 \xm 
PAH emission at 3.6 fim, which together contribute be- 
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tween 5-13% to the integrated light of the galaxies in 
our sample at 3.6 /xm and on average 20% locally to 
star-forming regions and 2) the spatial mapping of re- 
gions dominated by intermediate-age AGB or RSG stars, 
which contribute upwards of 50% to the 3.6 /xm flux at 
the locations of these isolated knots (but less than 5% to 
the integrated light of a galaxy). 

The hot dust and PAH at 3.6 /xm together are a po- 
tentially useful tracer of the thermal and non-thermal 
sources of dust heating in S 4 G galaxies (i.e. star for- 
mation and AGN), following calibration with ancillary 
8 /xm, 24 /xm, and/or Ha data available for part of the 
full sample. The contribution of intermediate-age stars 
at 3.6 /xm, which is separable from the dust emission 
via 8/xm flux (and possibly with uncorrected [3.6]-[4.5] 
color, as well), may also be useful for probing the diffuse 
component of 8 /xm non-stellar emission that does not 
trace actively star forming regions, but which is instead 
thought to originate with dust heating by an older pop- 
ulation of stars. In the future, we will use our maps of 
the regions dominated by intermediate-age AGB or RSG 
stars to constrain the fractional contribution of these ob- 
jects to emission in 3.6 /xm and blue- ward bands, as well 
as at 8 and 24 /xm. 

The limit of two input measurement bands enforced 
here, however, can challenge the detection of these 
sources with ICA, if only modestly. The intrinsically 
reddest regions, i.e. regions dominated by intermediate- 
age AGB and/or RSG stars, or HII regions where the 
flux in hot dust exceeds the spatially predominant PAH 
emission, in particular, tend to be over-subtracted. Since 
these are small, isolated zones that would be masked in 
modelling/structural decompositions in any case, we ar- 
gue that this is a small cost for the otherwise realistic 
correction throughout the rest of the contaminated area. 

In addition, we argue that errors associated with the 
ICA contaminant correction-i.e. the intrinsic incompati- 
bility between the number of real sources and the number 
of extracted sources-are not significant, the correction 
even in this case marking an improvement over alterna- 
tive corrections. Our use of ICA is superior to, e.g. the 
scaling and subtraction of the 4.5 /xm image from the 
3.6 /xm image, or to the subtraction of the non-stellar 
emission imaged at 8 /xm, since it can accommodate for 
the intrinsically different spectral shapes of the stars and 
dust between 3.6 /xm and either 4.5 or 8 /xm. 

In this way, our correction for contaminant emission at 
3.6 /xm and/or 4.5 /xm is optimal for a large data set such 
as S 4 G, which samples throughout the Hubble sequence, 
over a large range of dust contents, gas fractions, and 
star formation histories. Furthermore, the correction is 
self-reliant-using information from only the two IRAC 
bands available throughout the sample-and can assure 



uniformity in the treatment of contaminants for all ef- 
fected galaxies. Contaminants should be identifiable by 
deviations from the expected colors of K and M giants 
over significant areas in the original images, i.e. where 
colors appear in excess of [3.6]-[4.5]=0. 

Corrected 3.6 /xm images, which supply an optimal 
view of the underlying old stellar light, are ideal rela- 
tive tracers of stellar mass and should be indispensable 
for measuring structural diagnostics such as bar strength 
and ellipticity, spiral arm-interarm contrasts and rela- 
tive Fourier amplitudes and phases. Our maps can also 
be converted into absolute mass maps using any of the 
standard techniques for estimating the stellar M/L (e.g. 



Bell fc de Jong 2001 
optical and/or JN1R 



Zibetti et al.|2009 ) , where ancillary 
at a is available. 

will be explored in Paper II, 



On the other hand, as 
corrected images may support an independent view of 
the mass distribution in nearby galaxies. IRAC [3.6]- 
[4.5] colors, although more sensitive to photometric er- 
rors than those defined by a broader wavelength range, 
are virtually free of the systematic uncertainties plaguing 
the colors typically used to constrain the stellar M/L, e.g. 
due to dust extinction/reddening and young stars. In 
even the latest generation of m aps that actively co mpen- 



sate for these systematics (i.e. Zibetti et al. 



20091) there 



are lingering uncertainties: convergence on the TP- AGB 
phase of stellar evolution that dominates NIR light at ^1 
Gyr has yet to be reach ed among differe nt models (cf. 
Chariot fe Bruz ual 2007, |Maraston|[2005| and |Mouhcine 
|fc Langon||2002| ) ; where young and/or intermediate-age 
populations are identified, the M/L can be underesti- 
mated given an underlying old stellar population; al- 
though dust extinction is accounted for, the reflection 



al 



an d scat tering of stellar light by dust is not (Zibetti et 
_ J2009]). 
Li 



averaging our optimal view of the old stellar light we 
anticipate generating high-quality 2D mass maps that 
will provide an enhanced view of structure and mass 
in nearby S 4 G galaxies. Structural decompositions, 
measurements of bar and/or spiral arm torques, and 
improved constraints on the baryonic contribution to 
the distribution of mass in galaxies will be central to the 
insight on galaxy evolution that S G can deliver. 
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Observatory, which is a facility of the National Science 
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